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Abstract. We discuss the X-ray luminosity function and number counts of normal spiral galaxies. A detailed 
comparison is performed of the LFs compiled at infrared, radio and optical wavelengths and converted into XLFs 
using available relations in the literature with the XLF directly estimated in the 0.5-2 keV energy band from 
X-ray surveys (Norman et al. 2004). We find that the XLF from the local sample of IRAS galaxies (Takeuchi 
et al. 2003) provides a good representation of all available data samples; pure luminosity evolution of the form 
(1 + z) v , with r\ < 3, is favoured over pure density evolution. We also find that the local X-ray luminosity density 
is well defined, 3 ■ 10 37 ± 30% erg s" 1 Mpc -3 . 

We discuss different estimates of the galaxies Log iV-Log S, selected from the Chandra Deep Field surveys with 
different selection criteria, and find that these have similar slopes, but normalisations scattered within a factor 
~ 2, of the same order of the Poissonian error on the number counts. We then compare the observed galaxies 
Log TV-Log S with the counts predicted by the integration of our reference z = XLF. From the analysis of 
number counts alone, it is not possible to discriminate between density and luminosity evolution; however, the 
evolution of galaxies must be stopped in both cases at a redshift z ~ 1-2. The contribution of galaxies to the 
cosmic X-ray background is found to be in the range 6%-12%. By means of a complementary analysis with cosmic 
star formation models, we also find that the observed X-ray number counts might be not compatible with very 
large star formation rates at z ~ 3 as suggested by sub-mm observations in Blain et al. (1999a). 
Concerning the content of current and, possibly, future X-ray surveys, we determine the fraction of normal galaxies 
around the current flux limit. At ~ 3 • 10 -17 erg s _1 cm -2 in the 0.5-2.0 keV band, the normal galaxies in the 
Chandra Deep Field surveys are about the (30 ± 12)% of the total number of objects. At fainter fluxes the fraction 
of galaxies will probably rise, and overcome the counts from AGN at fluxes < 10 -17 erg s -1 cm~ 2 . 

Key words. X-rays: galaxies - galaxies: luminosity function - galaxies: evolution - galaxies: high-redshift - infrared: 
galaxies - galaxies: spiral 



1. Introduction 

The X-ray luminosity of normal spiral galaxies (hereafter 
'galaxies') appears to be a reliable, absorption-free estima- 
tor of star formation (i.e. their X-ray luminosity is dom- 
inated by emission from star formation related processes, 
and not by an AGN). Growing attention has been paid 
to its relation with emission at other wavelengths usu- 
ally taken as star formation indicators, such as the radio 
and far infrared (FIR) bands, in both the local and the 
distant universe; it has been found (Shapley et al. 2001; 
Bauer et al. 2002; Ranalli et al. 2003, hereafter RCS03) 
that the X-ray, radio and FIR luminosities are linearly and 
tightly correlated. Since the FIR and radio luminosities 
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are usually taken as Star Formation Rate (SFR) indica- 
tors (Kennicutt 1998; Condon 1992), it has been proposed 
that the X-ray luminosity also is a SFR indicator (RCS03). 

Thus, some authors (Nandra et al. 2002; Norman et al. 
2004) have begun to explore the X-ray emission as a tool 
to investigate the Cosmic Star Formation History (CSFH). 
To this end, the study of the X-ray luminosity function 
(XLF) of galaxies and of its evolution represents a neces- 
sary step. While the number counts may not lead directly 
to a determination of the CSFH, they are indeed useful 
as a testbed to constrain the evolution of the luminos- 
ity function (LF) at redshifts larger than those directly 
probed by the observations. 

Some early attempts to derive an XLF (Schmidt et al. 
1996; Georgantopoulos et al. 1999) were mainly limited 
by problems in selecting the galaxy samples. Nowadays, 
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the superb imaging capabilities of the ACIS detectors on- 
board Chandra has allowed the detection of extremely 
faint X-ray fluxes (<~ 10~ 17 erg s _1 cm~ 2 in the 0.5-2 
keV band), thus probing the X-ray emission from low lu- 
minosity AGN and galaxies at moderate redshifts (z < 1). 
However, in order to distinguish between accretion- and 
star formation-powered systems, spectroscopic and photo- 
metric follow-up observations are necessary. While these 
efforts are ongoing, the difficulty to unambiguously clas- 
sify the nature of the X-ray emission, makes a detailed 
determination difficult of the XLF properties rather hard. 

An attempt to overcome these difficulties, based on 
photometric selection criteria, has been discussed in 
Norman et al. (2004), where a sample was defined contain- 
ing 210 galaxies with known rcdshift from the Chandra 
Deep Field catalogues (Alexander et al. 2003; Giacconi 
et al. 2002). A Bayesian approach was chosen to derive 
a selection probability from the values of three different 
parameters: 

— the 0.5-2.0 keV X-ray luminosity for star forming 
galaxies is usually less than ~ 10 42 erg s _1 ; 

— star forming galaxies have a softer spectrum than 
AGN; this translates in selecting objects with an hard- 
ness ratio HR < -0.8 (with HR = (H - S)/(H + S), 
where H and S represent the 2.0-10 and 0.5-2.0 keV 
fluxes, respectively); 

— an X-ray/optical (R band) logarithmic flux ratio 
Log(Fx-ray/-Po P t) < —1 (a Log flux ratio of —1 is 
usually taken as an approximate boundary between 
galaxies and Seyferts; Maccacaro et al. 1988). 

Two redshift bins (z < 0.5 and 0.5 < z < 1.2 with 
mean redshifts z = 0.27 and 0.79, respectively) were con- 
sidered in order to have a comparable number of galaxies 
in both bins. A binned LF was derived with the method 
developed in Page & Carrera (2000), which is a variant of 
the classical 1/V max method by Schmidt (1968). The evo- 
lution of the XLF, parameterised with a single power law 
in both the redshift bins, is adequately described by a pure 
luminosity evolution of the form (1 + z) m with 771 ~ 2.7 
(error not available). It should be noted, however, that for 
a power law LF it is not possible to discriminate between 
density and luminosity evolution. 

A robust estimate of the XLF and of the cosmologi- 
cal evolution of star forming galaxies would greatly ben- 
efit from a detailed investigation of the relative merits 
of different selection criteria. For example, the Bayesian 
approach developed by Norman et al. (2004) might be af- 
fected by some contamination from low luminosity AGN 
with a soft X-ray spectrum. Well defined criteria to se- 
lect large samples of star forming galaxies are available in 
the infrared and radio bands, but the lack of X-ray ob- 
servations prevents a direct determination of their XLF. 
Nevertheless, it appears feasible to attempt an inter com- 
parison of the LFs obtained in the radio, mid/far infrared 
and optical bands for statistically significant samples, and 
in a rcdshift range as wide as possible, by transform- 
ing each of them in a corresponding XLF. This multi- 



wavelength approach builds on the well-known linear re- 
lations between X-ray, radio and far infrared luminosities 
calibrated for a well defined sample of objects, which we 
reported in RCS03, and also on non-linear relations be- 
tween blue and X-ray/radio/FIR luminosities reported by 
Shapley et al. (2001). 

Thus the structure of this article is as follows. In Sect. 2 
we review the current determinations of infrared, radio, 
and optical (B band) LFs. In Sect. 3, the LFs are converted 
into the X-ray band, and the predicted local (z = 0) XLFs 
are discussed; the predicted XLFs are also compared with 
the one observed by Norman et al. (2004). In Sect. 4 we 
review the observed X-ray Log A-Log S of galaxies based 
on different selection criteria. We recall from RCS03 the 
predicted X-ray counts based on the radio Log A-Log 5* 
and compare them with the observed X-ray counts. We 
also compare the observed counts with those predicted by 
the integration of the LFs. In Sect. 5, the Log A-Log S 
allow us to derive constraints on the LF evolution and on 
the cosmic star formation history. The contribution of the 
considered population of galaxies to the cosmic infrared 
background is also calculated as used to further constrain 
the LF evolution. Finally, in Sect. 7 we discuss the impli- 
cations for the XLF of star forming galaxies and for its 
evolution based on the previous analysis, we derive the 
fraction of galaxies in current deep X-ray survey and dis- 
cuss the implications for future surveys and missions. 

If not otherwise indicated (e.g., Table 1 where LF pa- 
rameters are given in an 7? - m dependent form), we assume 
H a = 70 km s" 1 Mpc" 1 , M = 0.3 and Q A = 0.7. 

2. The LFs of star forming galaxies 

We will consider the local differential luminosity func- 
tion ip(L)dLogL, which is the comoving number density of 
sources per logarithmic luminosity interval, subject to ei- 
ther a pure density evolution of the form ip(z) oc (l + z) r,d , 
or a pure luminosity evolution of the form L{z) oc (l+z) m . 

2.1. The FIR LFs 

Infrared surveys provide a powerful method to select star 
forming galaxies, since the bulk of the far and near infrared 
emission is due to reprocessed light from star formation, 
with AGN representing only a minor population (de Jong 
et al. 1984; Franccschini et al. 2001; Elbaz et al. 2002). In 
the following, we summarise the current determinations of 
the FIR luminosity functions. 

The FIR LFs discussed below can be safely assumed 
to be essentially unaffected by a contribution from Seyfert 
galaxies. By considering the 60/x LF of Seyfert galaxies 
(Rush et al. 1993), we estimate the fraction of Seyferts in 
the IRAS and ISO LFs discussed below to be ~ 5% for 
galaxies with vL v <~ 10 10 L©, and ~ 10% for galaxies with 
vL v ~ 1O 12 £0. 

Saunders et al. (1990) defined a sample of 2,818 galax- 
ies matching common selection criteria from different 
IRAS samples, with a flux limit around 0.6 Jy at 60fi, 
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Fig. 1. Infrared local luminosity functions. Solid curve: IRAS 
60/i LF from T03. Dotted curve and dark grey area: IRAS 
60/i LF from Saunders et al. (1990), with errors from 1/Vmax 
analysis. The two IRAS LFs coincide at bright luminosities. 
The size of the errors on the T03 LF and on the knee region 
of Saunders et al. is comparable with the line width. Dashed 
curve: ISO 90/t LF from S04; the shaded area corresponds to 
the la error bar. Dot-dashed curve: ISO 15/i LF from P04; the 
light grey area shows the la error bar. While the IRAS LFs are 
well constrained in the whole range of luminosities shown here, 
the ISO 90/i LF suffers from incompleteness for vL v < 1O 9 L0. 
The monochromatic luminosities at 90/i have been converted to 
bolometric ones by assuming vL v = cost. The 15/i luminosities 
have been converted by assuming L60m/^i5h ~ 5 (Mazzei et 
al. 2001). Ho = 70 is assumed. 

complete at the 98% level, redshifts z < 0.05, and de- 
rived a 60/i luminosity function. The best-fit model has 
the shape 



exp 
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-_Log 10 ll + - 



L 
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with the parameters reported in Table 1. The evolution 
was parameterised as pure density evolution and the ex- 
ponent found to be i]d = 6.7 ±2.3. The LF is shown as the 
dotted curve in Fig. 1, with errors from l/V majX analysis 
as the dark grey area. 

The 60/i luminosity function was revised by Takeuchi 
et al. (2003, 2004, hereafter T03) by enlarging the galaxy 
sample: 15,411 galaxies from the Point Source Catalog 
Redshift (PSCz, Saunders et al. 2000) were used, cover- 
ing 84% of the sky with a flux limit of 0.6 mJy at 60/i. 
The same parameterisation of Eq. (1) was used for the 
LF (Table 1). In Fig. 1 it is shown as the solid curve. 
This LF coincides with the Saunders et al. (1990) one 
for vL v > 10 9 L© and only differs by a factor of 2 at 
1O 7 L0. It was also found that the density evolution pro- 



posed by Saunders et al. (1990) is not consistent with this 
sample: a milder evolution, rj<i = 3.4 ± 0.70, is reported 
to be the best-fit description of the data. Since this rep- 
resents an update of Saunders et al. (1990) LF, in the 
following we shall refer only to the T03 LF. Note that 
the redshift distribution of the PSCz galaxies (which is a 
superset of Saunders et al. sample) spans a very limited 
range (z < 0.07), so that the evolution exponent is poorly 
constrained. 

The LF at 90/i was recently determined by the ELAIS 
survey team (Serjeant et al. 2001, 2004, hereafter S04) 
with the PHOT instrument onboard the ISO satellite: 151 
objects detected in an area of 7.4 square degrees with flux 
limit 70 mJy and spectroscopic redshifts (most of the ob- 
jects having z < 0.3) were used to derive the LF. It has 
not been fitted to a parametric form; see the dashed curve 
in Fig. (1) for the LF and Table 1 for the evolution. To 
convert the ELAIS LF to 60/i, a Seofi/Sgo^ — 0.66 flux 
ratio can be derived by assuming vh v = cost, which well 
reproduces the spectrum of star forming galaxies in the 
wavelength region here considered. 

2.2. The mid infrared LF 

We consider the recent determination of the mid-infrared 
(MIR) LF by Pozzi et al. (2004, hereafter P04) with 
ISOCAM data at 15/i from the ELAIS survey. The local 
LF was derived from a sample of 150 objects with z < 0.4 
and flux <Si5 M > 1 mJy observed in two fields for a total 
area of 2.1 square degrees. P04 divided the galaxies be- 
tween quiescent spirals and starburst on the basis of the 
MIR to optical flux ratio. The best-fit LF is given as two 
Saunders functions (Eq. 1) with the parameters shown in 
Table 1. We converted these LFs to the 60/i wavelength 
using the observed relation Lqq^/ ~ 5 from Mazzei 
et al. (2001). The sum of the two 60/i LFs is shown in 
Fig. (1) as the dot-dashed line. The grey area shows the 
la error bar derived from a 1/V max analysis. This local LF 
matches well the IRAS ones in the region around the knee 
(i.e. 10 10 < L < 10 11 LQ), while at lower and higher lumi- 
nosities it stays a factor of two-three lower (see discussion 
in P04); the evolution parameters were determined both 
by analysing the LF in the redshift interval in which it is 
defined and by fitting the MIR source counts. P04 found 
that the evolution of the MIR LF can be described as fol- 
lows: no evolution for the spirals from z = to ~ 0.4, and 
both luminosity and density evolution for the starbursts 
with z < 1 (no evolution for the starbursts beyond redshift 
1). The best fit parameters for evolution are rji 
and ?7d = 3.8 ± 2. 



2.3. The radio LFs 

Among the first determinations of the radio (1.4 GHz) 
LF of galaxies, the LF by Condon (1989, hereafter C89) 
was derived by considering the radio observations of 292 
spiral galaxies with blue magnitude B < 12 and declina- 
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Table 1. Luminosity functions of galaxies at infrared, radio and blue wavelengths. 



rJanu 


Author 


LF form 


Parameters 


Evolution 




Saunders 
et al. 


Eq. (1) 


<p* = (2.6 ± 0.8) • 10~ 2 h 3 Mpc -3 , a = 1.09 ± 0.12, 
cr = 0.72 ± 0.03 and L* = (2.95t?;!£) • 10 8 h~ 2 LQ 


r; d = 6.7 ±2.3 


60ri 


Takeuchi 
et al. (T03) 


Eq. (1) 


<p* = (2.34 ± 0.30) ■ 10~ 2 h 3 Mpc~ 3 , a = 1.23 ± 0.04, 
a = 0.724 ± 0.01 and L* = (4.4 ± 0.9) • 10 8 h~ 2 LQ 


r? d = 3.4 ± 0.7 


90/i 


Serjeant 
et al. (S04) 


numerical, see Fig. 


(1) 


r? d = 3.5 ±1.1 


15fi 


Pozzi 

et al. (P04) 


Eq. (1) spirals ± 
Eq. (1) starbursts 


quiescent spirals: 

<p* = 3.55 ■ 10~ 3 h 3 Mpc~ 3 , a = 1.10 ± 0.25, 
a = 0.5±g;| and L* = (6.3±j? 5 5 ) • 10 s h~ 2 LQ 
starbursts: 

ip = 2.95 • 10 h Mpc , a fixed at 0, 

o = 0.39 ± 0.025 and L* = (6.3±|;l) • 10 8 h~ 2 LQ 


spirals: no evolution 

starbursts: 
ri d = 3.8 ± 2 
with 

m = 3.5tj : g 


1.4GHz 


Machalski & 

Godlowski 

(MGOO) 


Eq. (1) 


ip* = (7.9t 3 ; 8 ) • 10~ 3 h 3 Mpc^dex" 1 , 

a = 1.22 ± 0.27, a = 0.61 ± 0.05 and 

Log L* = (28.13 ± 0.37) h~ 2 erg s" 1 Hz" 1 


constraints are too loose 
to determine the evolution 


1.4GHz 


Condon (C89) 
+ Haarsma 
et al. (H00) 


Eq. (2) 


Y = 2.85, X = 22.40, W = 2/3, B = 1.5 
(assuming Ho = 50) 


lfi(z) CX (1±2)" ' 0579 

• ex p[-(Tfs) 231 ] 
with 

L(z) oc (1±2) 3 ' 97 

•^[-(Tig) 1 ' 02 ] 


1.4GHz 


Sadler, 

et al. (S02) 


Eq. (1) 


<p* = (7.8 ± 0.7) ■ lO" 2 h 3 Mpc- 3 dex- 1 , 
a = 0.84 ± 0.02, a = 0.940 ± 0.004 and 
Log L* = (27.15 ± 0.03) h~ 2 erg s" 1 Hz" 1 


constraints are too loose 
to determine the evolution 


Blue 


Wolf et al. (W03) 


15 Schechter functions with different parameters according to spectral type and redshift bin 



tion north of 6 = —45° from the Revised Shapley-Ames 
(RSA) catalogue (Sandage & Tammann 1981). The C89 
LF, shown in Fig. (2) as the solid line, is parameterised as 

Log</?(LogL) = 28.83 + Y- 1.5 LogL 

_^ + (^x£y' 2 Mpc -3 dS3t -, (2) 

with the best-fit parameters shown in Table 1. In addition 
C89 considered also a sample of 142 galaxies selected from 
the Uppsala General Catalogue of Galaxies (UGC, Nilson 
1973) with declination —2.5° < 5 < 82°, angular diameter 
> 1.0' and classified as "starburst" or AGN on the basis of 
their radio morphology, infrared/radio and 60^/25^ flux 
ratios. It was found that the LF parameters derived for 
the RSA spirals fitted well also the UGC starbursts. 

The most recent determinations of the radio LF of 
star forming galaxies are based on the cross correlation 
of galaxy redshift surveys at optical wavelengths with ra- 
dio all-sky surveys. First we consider the sample of radio 
detected galaxies in the Las Campanas Redshift Survey 
(Machalski & Condon 1999), which comprises 1157 galax- 
ies with magnitude R < 18.0 detected in the NRAO VLA 
Sky Survey (NVSS), i.e. brighter than 2.5 mJy beam -1 



at 1.4 GHz: 502 galaxies were classified as "starburst" 
and 655 as "AGN" on the basis of their FIR-radio flux 
ratios, 25/x-60^i spectral indexes and radio-optical flux ra- 
tios. Machalski & Godlowski (2000, hereafter MGOO) de- 
rived the LFs for both classes of objects; here we consider 
only the "starbursts" (Table 1; dashed line in Fig. 2). 

The redshifts of the Las Campanas "starburst" galax- 
ies are in the range 0.01 < z < 0.075. As in the case of the 
IRAS galaxies, this redshift range is too small to allow a 
reliable estimate of the evolution. The constraints on the 
evolution found by MGOO are too loose to either confirm 
or reject cosmic evolution for the Las Campanas galaxies. 

We also consider the radio LF of galaxies from the 2dF 
Galaxy Redshift Survey (2dFGRS; Sadler et al. 2002, here- 
after S02): 242 galaxies with redshift z < 0.3 were identi- 
fied in the cross-correlation of 2dFGRS sources (covering 
an area of 325 deg~ 2 ) with the NVSS. The galaxies were 
classified as star forming systems (different from AGN) on 
the basis of their optical spectra. To further extend the LF 
at luminosities lower than ~ 10 29 erg s _1 Hz -1 , S02 in- 
cluded also the RSA spirals previously considered by C89. 
The LF was parameterised as a Saunders function, with 
the best fit parameters shown in Table 1; it is the dot- 
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Fig. 2. Radio local luminosity functions. Solid curve: LF from 
spiral galaxies in the RSA (C89); light grey area: error on 
the RSA spirals LF from 1/Vmax analysis. Dot-dashed line: 
radio LF from star forming galaxies in both the 2dF Galaxy 
Redshift Survey and the RSA spirals (S02); dark grey area: er- 
rors from 1/V m ax for the 2dFGRS galaxies. Dashed curve and 
diagonally shaded area: LF for starburst galaxies in the Las 
Campanas Redshift Survey (MGOO) and corresponding 1 / V ma x 
error. Ho = 70 is assumed. 



dashed line in Fig. (2). A V/V max test for evolution of 
the 2dFGRS galaxies yielded inconclusive results, mainly 
because of poor statistics (at the flux limit of the NVSS, 
only the highest luminosity bins of the LF are sampled for 
0.2<z<0.3). 

The behaviour of the radio LF at higher redshifts 
(z < 2) was explored by Haarsma et al. (2000), who built 
a sample of galaxies detected in deep radio surveys. They 
considered the C89 local LF and fitted its evolution to 
the high redshift data. Haarsma et al. (2000) considered 
a rather complex form for the evolution with six free pa- 
rameters (see Table 1). However, their best-fit evolution 
parameters closely resemble a pure luminosity evolution 
not too far from the one found by T03, with the added 
constraint of an exponential cutoff around z ~ 1.4. In the 
following, we shall refer with 'C89 LF' to the C89 LF with 
the evolution determined in Haarsma et al. (2000). 

2.4. The blue LF 

The evolutions of the infrared and radio LFs discussed 
in the previous section are quite loosely constrained, ei- 
ther because of the limited redshift range, or because 
of the limited number of objects present in the surveys. 
However, larger surveys are currently ongoing at optical 
wavelengths exploring wider redshift ranges. Here we con- 



sider the COMBO-17 survey (Wolf et al. 2003, hereafter 
W03) which provides a blue (Johnson B-band) LF from 
25,000 galaxies with 0.2 < z < 1.2. Redshifts are obtained 
from photometric observations taken with a set of 17 fil- 
ters. The foremost data analysis goal of the COMBO-17 
approach is to convert the photometric observations into 
a low resolution spectrum that allows simultaneously a re- 
liable spectral classification of stars, galaxies of different 
types and quasars, as well as an accurate redshift estimate. 

The galaxy B-band LFs are given in W03 for five red- 
shift bins (from z = 0.2 to 1.2), so that no assumptions 
have to be made about evolution. The LFs are also given 
for four different spectral types defined on the basis of the 
grid of template spectra by Kinney et al. (1996). Since 
spectral types may be broadly matched to morphological 
types (Madgwick et al. 2002), we consider the LFs (one 
for each redshift bin) for late-type galaxies as derived by 
W03, i.e. with the exclusion of spectral type 1 in their 
classification (see Table 1). 

3. The predicted XLFs at z = and their 
evolution 

The LFs discussed in Sect. 2 are converted by using the 
approach first developed in Avni & Tananbaum (1986) 
(see also Georgantopoulos et al. 1999 and Norman et al. 
2004): 

+ 00 

^x(Log L x ) = J ^v(Log L Y ) 

-co 

■P (Log L x |Log L Y ) dLog L Y (3) 

where P (Log Lx|Log Ly ) is the probability distribution 
for observing Lx for a given luminosity in the Y band. In 
RCS03 it was reported that the X-ray luminosity is tightly 
correlated with radio and FIR luminosities (LogL .5-2 = 
LogLi.4GHz + 11.08 and LogLn.5-2 = LogL 60M + 9.05, with 
£0.5-2 in erg s _1 , and with L1.4 and L 60fJ , in erg s _1 Hz^ 1 ). 
By assuming a Gaussian probability distribution for these 
correlations, one has 



P(Log £0.5-2! Log L 60ai ) = 



Log I/ 6 Q M +9.05-Log I/p.5-2 



27TO- 



(4) 



with a ~ 0.30, and for the radio band 

P(Log L .5-2|Log £i.4GHz) = 

I Log L 1 4 + II.IO — Log Lq 5 _ 



2ir<r 



(5) 



with a ~ 0.24. 

To convert the W03 LF into an XLF we consider the 
B-band/X-ray correlation, which has been the subject of 
much attention in X-ray studies of galaxies of all morpho- 
logical types (Giacconi & Zamorani 1987; Fabbiano 1989; 
David et al. 1992; Shapley et al. 2001; Fabbiano & Shapley 

2002). A correlation with the form Lx oc L B 



+0.1 
2 was first 
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found by Fabbiano et al. (1988) by analysing 51 spiral 
galaxies observed with Einstein. An enlarged sample of 
234 spirals from the Einstein archive was built by Shapley 
et al. (2001), from which we derive L x cx L^ 13±0 09 by 
performing a linear regression analysis (Isobe et al. 1990) 
with the Schmitt method in the ASURV code for censored 
data analysis (La Valley et al. 1992), which we chose be- 
cause of the large number of upper limits on the X-ray 
luminosities and consistently with Shapley et al. (2001). 
Thus we use the following relationship 



P (Log Lo.5 



- 2 |Log L B ) 
1 



1.13Log L B -8.97-Log L B _ 2 



27TO- 



(6) 



with a ~ 0.54. 

We note that while the B/X-ray correlation was de- 
rived by making use of magnitudes corrected for extinc- 
tion, caused by the inclination of the galaxy on the plane 
of the sky, the W03 LF is not corrected for extinction. 
However, this should not be a serious problem, since 
the median inclination-induced extinction for the galaxies 
in the Third Reference Catalogue of Galaxies (RC3, de 
Vaucouleurs et al. 1991) is 0.26 magnitudes (i.e. a factor 
of ~ 1.3), so that we have to correct the LF for inclination- 
induced extinction by rescaling the normalisation of the 
LF by a factor 1.3. 

The radio and IR derived XLFs are shown in Fig. (3). 
Around the 'knee' (i.e. 10 40 < L x < 10 41 ) all the predicted 
XLFs agree within a factor of 2, which is the approximate 
size of the error bars on the radio and infrared LFs. Much 
more scatter is present among the tails of the LFs, due to 
the small number statistics. 

The main point here is that the infrared and radio 
LFs define a locus where an observed local XLF should 
fall. This confirms the effectiveness of the various selec- 
tion criteria of galaxies used in the surveys. The X-ray 
luminosity density (XLD) of galaxies in the local universe 
(XLD = f_°° dLogL ip(LogL) ■ L) is also well defined, with 
estimated values that are spread within a factor ~ 1.7 
around a formal average ~ 3.4 • 10 37 erg s _1 Mpc -3 ; the 
T03 XLD being close to the average value (Table 2). Thus 
the T03 XLF, built on the largest sample of galaxies, ap- 
pears well suited to be representative of the local XLF 
(moreover, had we weighted the mean by using the num- 
ber of galaxies in surveys, the average XLD would be co- 
incident with the T03 one). 

We also computed the percentage contribution to the 
cosmic X-ray background (CXB) in the 0.5-2.0 keV en- 
ergy band by adopting the mean XLD and a uniform dis- 
tribution (no evolution) of the galaxies up to z — 4; as 
regards the CXB we have adopted the spectral shape of 
9.8_E~ 14 photon s _1 cm~ 2 sr _1 keV -1 , with normalisa- 
tion taken from Revnivtsev et al. (2005). By adopting the 
evolution parameters as discussed in Sect. 4 we find that 
the contribution to the CXB is comprised in the percent- 
age interval 6%-12%, keeping in mind, however, that the 
average spectral slope of the galaxies is much softer than 
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Fig. 3. (In colour only in the electronic edition) IRAS, ISO, 
and radio local luminosity functions converted to the X-rays. 
Solid (black) curve: IRAS 60/i LF from T03. Continuous line 
with triangles: ISO QQa LF (S04). Dot-dashed (blue) curve: 
ISO 15 fi LF (P04). Long-dashed (red) curve: radio LF (MGOO). 
Short-dashed (dark red) curve: radio LF (C89). Short-dash- 
long-dash (green) curve: radio LF (S02). 

Table 2. X-ray luminosity densities of galaxies at z = in the 
0.5-2.0 keV band. 



Luminosity function 



Predicted luminosity density 
(erg s -1 Mpc -3 ) 



Takeuchi et al. (T03) 


3.0 


10 37 


Condon (C89) 


4.1 


10 37 


Sadler et al. (S02) 


3.7 


10 37 


Pozzi et al. (P04) 


2.5 


10 37 


Serjeant et al. (S04) 


4.1 


10 37 


Machalski & Godlowski (MGOO) 


2.9 


10 37 


average value 


3.4 


10 37 



that of the adopted extrapolation of the CXB. These find- 
ings are consistent with the 1% contribution to the XRB 
by galaxies detected in the CDFs (Hornschemeier et al. 
2003). 

3.1. The evolution of the XLFs 

The infrared-, radio- and blue-derived LFs are compared 
with the X-ray observed LF in Figs. (4) and (5), where all 
the LF are drawn only in the X-ray luminosity range cor- 
responding to that spanned by the data defining each LF 
in the original band. Since the Norman et al. (2004) XLF 
was derived at higher redshifts (z ~ 0.27 and z ~ 0.79), 
in order to perform a meaningful comparison we have to 
consider the different constraints and parameterisations 
on the evolution of the LFs. The LFs whose evolutions 
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Fig. 4. Comparison of the X-ray luminosity functions, (a) Low-z bin [z = 0.27), density evolution, (b) Low-z bin (z = 0.27), 
luminosity evolution, (c) High-z bin [z = 0.79), density evolution, (d) High-z bin (z = 0.79), luminosity evolution. Filled squares 
or circles: observed XLF from Norman et al. (2004). Solid curve: IRAS LF (Takeuchi et al. 2003, 2004); the grey shading shows 
the la error on the evolution exponent; triangles within line shading: ISO LF (Serjeant et al. 2004). The XLFs have been evolved 
with the prescriptions listed in Table 1 to the median redshifts of the two bins in which the observed XLFs were obtained. 



are constrained only by local data [z < 0.3, i.e. the T03 
and S04 LFs) are shown in Fig. (4). Only pure density 
evolution was considered for these LFs by their respective 
authors. Nonetheless, it seems that pure luminosity evolu- 
tion may provide an equally good fit to the data, at least 
for the T03 LF for which the analysis is ongoing (Takeuchi, 
priv. comm.). However, we have decided to explore also 



the effect of pure luminosity evolution by adopting the 
same exponents reported for the density evolution (this 
is not formally correct but may serve as a first-order ap- 
proximation). The other LFs (P03, H00 and W03), whose 
evolutions are constrained from higher redshift data, are 
shown in Fig. (5). These are either parameterised in both 
luminosity and density (P04 and H00), or are observed 
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Fig. 5. Comparison of the X-ray luminosity functions. Left panel: low-z bin (z — 0.27). Right panel: high-z bin (z = 0.79). Data 
points: observed XLFs from Norman et al. (2004). Solid curve: ISO XLF (P04), with la error on evolution (see text); dashed 
curve: radio XLF (C89); dot-dashed curve: blue XLF (W03), observed in the redshift bins 0.2-0.4 and 0.6-0.8, respectively. The 
MIR and radio XLFs have been evolved with the prescriptions listed in Table 1 to the median redshifts of the two bins in which 
the observed XLFs were calculated. 



LFs in given redshift bins (W03). Two LFs arc not fur- 
ther discussed: Saunders et al. (1990), because it is indis- 
tinguishable from the T03, and the S02 and MG00 LFs 
because the constraints on their evolution are too loose to 
add significant information to our discussion. 

In both figures the shaded areas show, where available, 
the la error on evolution, while the errors on the z = 
LF are ignored. 

Considering pure density evolution (Figs. 4a and 4c), 
we find that the X-ray derived LF is in agreement with 
the FIR- and radio-derived XLFs for L x < 10 415 , but 
it is significantly flatter at its bright end (Lx > 10 415 ), 
where it has a density of galaxies which is about one or- 
der of magnitude larger than that of the FIR-radio galax- 
ies. This discrepancy is partly removed when we consider 
pure luminosity evolution (Figs. 4b and 4d), as it remains 
only noticeable (panel b) for low-redshift, high-luminosity 
(Lx > 10 42 erg s _1 ) galaxies. Although this discrepancy 
is only within 2a, it may be symptomatic of a residual 
fraction of low luminosity AGN embedded in the galaxies 
at the brightest end of the X-ray sample used by Norman 
et al. (2004). 

The XLF derived from the MIR LF (P04) appears to 
lie below the observed XLF in the low redshift bin (Fig. 5, 
left panel) and to cross the XLF in the high redshift bin 
(same figure, right panel). This might be partly explained 
because of the decoupling between the evolution of quies- 
cent and actively star forming galaxies: the net effect is 
that only the high luminosity (Lx > 10 40 erg s _1 ) tail 



of the LF evolves, but the non-evolving, lower luminosity 
end of the LF cannot fit the lowest luminosity, low redshift 
point of the observed XLF. 

The XLF derived from the radio LF (H00) crosses 
within la most of the observed XLF data points, with 
the exception of the two highest luminosity points in the 
low redshift bin (Fig. 5, right panel). 

The XLF derived from the blue LF (W03; dot-dashed 
curve in Fig. 5) gives results which are very similar to the 
T03 XLF. 

Thus, it seems that pure luminosity evolution is some- 
what favoured over pure density, with r/i between 2.7 
(the best fit value found by Norman et al. 2004) and 3.4 
(adopted from T03). 

The W03 XLF offers the possibility to check the evo- 
lution independently from the Norman et al. (2004) XLF 
and in a larger number of redshift bins. Moreover, be- 
cause of the large size of the galaxy sample and the wider 
redshift range, the luminosity and density evolution can 
be analysed separately. W03 reports on the evolutionary 
patterns for each spectral type: on the whole, these pat- 
terns can be understood as an increase in starburst activ- 
ity with redshift, with an activity maximum at z > 1, in 
combination with a propagation of the galaxies through 
the spectral types as the mean age of the stellar popula- 
tion decreases. Quantitatively, we may consider the red- 
shift evolution of both the X-ray luminosity density and 
of the number density of the W03 galaxies. We find that 
the luminosity density doubles from z ~ 0.3 to 1.1, while 
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Fig. 6. Comparison of the T03 and W03 XLFs for each red- 
shift bin in which the W03 LFs are denned. Pure luminosity 
evolution with r\\ — 3.4 ± 0.7 is assumed for the T03 XLF. 
Solid line with grey error area: T03, with la error on evolu- 
tion. Dotted line: W03. The XLFs are progressively shifted for 
clarity by 1 dex along the vertical axis. The agreement is quite 
good in the low redshift bins and worsens at higher z. A better 
agreement is found if one assumes 771 = 2.7 for the T03 XLF, 
corresponding to the lower bound of the grey error area. 



the number density stays almost constant. This pattern is 
consistent with luminosity evolution prevailing over den- 
sity. 

A detailed comparison of the W03 XLFs with the 
evolved T03 XLF (pure luminosity, r\\ = 3.4) shows a fairly 
good agreement (less than a factor of 2) in the redshift bins 
from z ~ 0.3 to ~ 0.7, while it is a factor of ~ 3 below the 
T03 XLF in the z ~ 1.1 bin (Fig. 6). No agreement would 
occur had we considered pure density evolution. The as- 
sumption of pure luminosity evolution for the T03 XLF, 
which was arbitrarily introduced in Sect. 3, appears to be 
validated by this comparison. We notice, however, that 
the evolved T03 systematically exceeds the W03 XLF for 
z > 0.7, perhaps an indication that the 771 = 3.4 evolu- 
tion is too strong; in fact, a milder exponent (e.g. 771 = 3) 
would reconcile the T03 and W03 XLFs for z < 0.9. In 
the z ~ 1.1 redshift bin the T03 XLF would still exceed 
the W03 one by a factor of 2-3, possibly hinting that the 
evolution should be stopped at z ~ 1. Altogether the two 
by far largest samples of sources provide a complemen- 
tary and self consistent picture up to a redshift z ~ 1, and 



in the analysis performed in the following Sects, we shall 
only consider the T03 XLF. 

4. The X-ray Log A^-Log S of galaxies 

We now turn our attention to the analysis of the X-ray 
number counts of galaxies, since they represent an imme- 
diate and alternate test bed for both the selection criteria 
in X-ray surveys and for the shape and evolution of dif- 
ferent LFs. 

4.1. The observed Log N-Log S 

The normal galaxy number counts at faint fluxes have 
been presented several times in the literature. While we 
just mention the first prediction based on optical counts 
by Giacconi & Zamorani (1987), we focus here on the dif- 
ferent determinations of the galaxy Log iV-Log S based 
on Chandra Deep Fields catalogues, which are defined at 
fluxes 3 • 10~ 17 < So.5-2keV < 1CT 15 erg s" 1 cm" 2 . 

The first determinations were made by Bauer et al. 
(2002), a sample of 11 galaxies selected by cross correlating 
X-ray and radio sources with optical spectra classified as 
"emission line galaxies" in the inner 3' of the CDFN, and 
by Hornschemeier et al. (2003), a sample of 43 galaxies 
with Log(Fx-ray/-Fopt) < —2.3 (a more stringent criterion 
to avoid AGN contamination). Hornschemeier et al. (2003) 
found also that the fraction of early type systems in their 
sample is somewhat low (~ 15%); thus we believe that 
in the following analysis the contribution form early type 
galaxies may be safely ignored. 

By refining the selection criteria it has been possible 
to construct the Log A^-Log S for larger samples to which 
we shall refer in our discussion. These are: 

— The Log A^-Log S of the Bayesian sample defined in 
Norman et al. (2004) and discussed in Sect. 1. By con- 
sidering the sky coverage of the Chandra Deep Field 
North (CDFN) survey, as reported in Alexander et al. 
(2003), we derived the Log A^-Log S for these sources, 
which is shown in Fig. (7) as the long-dashed his- 
togram; 

— Two samples defined in Bauer et al. (2004). The first 
one comprises 175 sources with internal absorption 
Ah < 10 22 cm -2 and hardness ratio HR< —0.8 and 
X-ray luminosity Lo.5-8keV < 3 • 10 42 erg s , or with 
off nuclear emission 1 . The second one comprises 109 
sources with Log(Fx- r ay/^ 1 o P t) < —1 and < 10 22 
cm~ 2 and HR< -0.8 and L .5_ 8 keV < 3 • 10 42 erg s" 1 
and no broad lines in the optical spectrum. Because 
of their different sizes, the two samples are named 
as 'optimistic' and 'pessimistic' estimates, respectively 
(Fig. 7). These samples also represent an update of the 
Hornschemeier et al. (2003) analysis (Hornschemeier, 
priv. comm.). 

1 A source is catalogued as 'off nuclear' if point-like X-ray 
emission is found within the optical contour of a galaxy in an 
off-centre position. 
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Fig. 7. (In colour only in the electronic 
edition) Observed X-ray counts. The 
thick, upper line and the horn-shaped 
symbol show the observed Log iV-Log S 
for all X-ray sources in the Chandra 
Deep Fields (Moretti et al. 2003) and 
the limits from the fluctuation analy- 
sis (Miyaji & Griffiths 2002a,b). Short- 
dashed red histogram and grey area: 
sources with Log (X-ray/optical flux ra- 
tio) < —1 with Poissonian error (this 
work). The width of the error area may 
be taken as representative for the fol- 
lowing Log iV-Log S too. Long-dashed 
black (upper) histogram: sources from 
the Bayesian sample of Norman et al. 
(2004). Continuous blue histograms: 
'optimistic' (upper) and 'pessimistic' 
(lower) estimates from Bauer et al. 
(2004). The data points and the dot- 
dashed green histogram at fluxes > 
10~ 15 erg s _1 cm -2 show the Log N- 
Log S of serendipitous spiral galaxies in 
XMM-Newton (Georgakakis et al. 2004) 
and ROSAT HRI fields (Tajer et al. 
2005), respectively. 



— A sample of 142 sources from the CDFN and the 
Chandra Deep Field South, which we selected for hav- 
ing an X-ray/optical flux ratio < — 1. Although this se- 
lection criterion may be regarded as somewhat crude, 
it is indeed useful since it allows an immediate ob- 
ject selection with the least possible number of param- 
eters. The resulting Log A-Log S is plotted as the 
short-dashed histogram in Fig. (7). An estimate of the 
Poissonian error (after Gehrels 1986, with confidence 
level of 68.3%) for this Log A^-Log S is also shown as 
the grey area. The width of the error area may also be 
taken as representative for the other Log A-Log S. 

We believe that our choice, although not exhaustive, 
characterises both the scatter among different number 
counts determinations and the 'state of the art' selection 
criteria. We find that there is a quite good agreement be- 
tween the slopes of the observed counts, while the normal- 
isations are comprised within a factor ~ 2. 

These results are somewhat strengthened by the re- 
cently determined galaxy counts at brighter X-ray fluxes 
(£0.5-2 ^ 10~ 15 ). A smooth transition is observed be- 
tween the bright tail of the number counts of the galaxies 
in the deep fields and the recently determined Log N- 
Log S of galaxies with bright X-ray fluxes, as shown 
in Fig. 7: Georgakakis et al. (2004) selected 11 sources 
from the cross correlation of serendipitous X-ray sources 
in archival XMM-Newton fields with galaxies from the 



Sloan Digital Sky Survey (SDSS); while Tajer et al. (2005) 
found 32 sources (from which we selected 21 galaxies 
with late-type morphology) by cross correlating a sample 
of serendipitous X-ray sources in archival ROSAT HRI 
fields with galaxies in the Lyon-Meudon Extragalactic 
DAtabase (LED A). The observed integral counts appear 
to be consistent with one power law (exponent ~ —1.4) 
spanning about four decades in flux. For easy reference 
in Figs. (8- 9) we always reproduce the counts obtained 
in the Chandra Deep Fields, but omit the counts at the 
brighter hand that would require unnecessary larger fig- 
ures' size. 

4.2. X-ray counts from integration of the LFs 

The observed X-ray counts may be checked against the 
number counts obtained by integrating the luminosity 
functions previously derived. In the following, we shall 
only discuss the number counts obtained from the T03 
XLF, since it has already been shown that its properties 
are the best defined ones and consistent with all the XLFs 
derived from the other samples. 

The number counts may be obtained by: 

Z„ax imax 

A(> S)= J dz J dLogL tp(LogL, z) ^ (7) 
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where dV/ dz is the comoving volume between z and z+dz. 
The K-correction can be safely neglected because the 
mean slope of the X-ray photon spectra of star forming 
galaxies is T ~ 2.1 (RCS03). We perform the integration 
in the luminosity interval 10 39 < Lx < 10 43 erg s" 1 at 
z = 0; when considering luminosity evolution, the inte- 
gration limits are also 'evolved', i.e. they are shifted by 
(1 + z) m . We note that, had we not evolved the integra- 
tion limits, the main conclusions would be unchanged; the 
only noticeable difference would be a larger number (by 
a factor <~ 3) of very faint sources (flux < 10~ 18 erg s _1 
cm -2 ). We perform the integration in the < z < 2 red- 
shift range by assuming two different possibilities for the 
evolution: the full range up to z < 2 and a restricted range 
up to z < 1. 

In Fig. (8, left panel) we plot the counts obtained from 
the T03 XLF with density evolution (r) d = 3.4). The in- 
tegrated counts lie on the lower side of the Log N-Log S 
bundle and are consistent with the 'pessimistic' estimate 
in Bauer et al. (2004). 

On the other hand, if we assume pure luminosity evo- 
lution with 771 = 2.7 as determined in Sect. 3.1 (Fig. 8, 
right panel), the integrated counts would be consistent 
with the observed ones, staying in the middle of the ob- 
served counts bundle for So. 5-2 > 5 • 10~ 17 erg s" 1 cm~ 2 . 
Performing the luminosity evolution up to z ~ 1 or 2 
would result in a difference of less than a factor of 2 at 
fainter fluxes, with the predicted counts at 3 • 10~ 17 erg 
s _1 cm~ 2 falling in the middle or in the upper end of the 
bundle, respectively. Had we assumed a stronger exponent 
{rj\ = 3.4 as explored in Sect. 3.1), the predicted counts 
would have outnumbered the observed ones. 

As a final comment we notice that the scatter in the 
observed Log TV-Log S (both the systematic scatter be- 
tween counts derived with different selection criteria and 
the Poissonian error on the single Log iV-Log S, Fig. 7) 
does not allow in itself to disentangle between density or 
luminosity evolution, but the above results strengthen the 
indication in favour of a luminosity evolution with rj\ <~ 2.7 
already emerged in the discussion of the XLFs of Sect. 3.1. 

We have so far explored the expected behaviour of the 
Log iV-Log S by integrating up to z = 2. However, it is 
known that star forming galaxies do exist and emit X-rays 
at redshifts z ~ 3, at least in the form of Lyman-break 
galaxies (Nandra et al. 2002). Nonetheless, our parame- 
terisation may still be justified as very little is known of 
the cosmic evolution of galaxies beyond z ~ 1-2, and it 
is not excluded that galaxies become less luminous, which 
for practical purposes would not be too far from the case 
hereby considered. 

In order to investigate the number counts from galax- 
ies at larger redshifts we consider the following scenario: 
galaxies evolve in luminosity up to a redshift z cv , and 
evolve no more beyond, while the integral is performed 
up to z = 4. By assuming r/i = 2.7, we find that the in- 
tegrated counts are consistent with the observed ones if 



A similar argument could be made for pure density 
evolution: in this case, the observed all-inclusive X-ray 
counts would not be outnumbered, but in order not to 
exceed the cosmic infrared background, one must have 
z cv < 2 (see Sect. 5). 

5. X-ray counts and the cosmic star formation 
history 

A first attempt in using the evolution of the XLF to put 
constraints on the cosmic star formation history (CSFH) 
was made in Norman et al. (2004), where the observed 
X-ray luminosity density was converted to a SFR density 
at the redshifts z = 0.27 and 0.79. Here we investigate 
a different possibility, which is to use the X-ray number 
counts to constrain the CSFH at higher redshifts (z < 4). 

The shape of the CSFH beyond z ~ 1-2 is still de- 
bated. Results from UV observations (Lilly et al. 1996; 
Madau et al. 1996) pointed towards a peak around z ~ 2 
with a SFR(z = 2) ~ 10 times that at z = 0. Pettini et 
al. (1998) applied a correction for dust extinction, which 
raised the peak SFR by a factor 2 and at the same time 
shifted the peak redshift at z <~ 3. However, sub-mm 
observations with the SCUBA camera could not be rec- 
onciled with the UV-derived, extinction corrected CSFH, 
leading Blain et al. (1999a) to propose the existence of a 
new population of luminous, strongly obscured, high red- 
shift (z < 5) galaxies. The CSFH derived from sub-mm 
observations have a peak at z ~ 3 with a SFR <~ 100 times 
the current (z = 0) one; at higher redshifts the data allow 
the SFR to stay constant or decline. 

In Sect. 4.2 we have discussed a particular case, in 
which the galaxies are distributed according to the evolved 
T03 XLF in the redshift interval 0-4, and we found that 
the XLF is not allowed to evolve beyond z ~ 1-2 in or- 
der not to violate the observed X-ray source counts. This 
would be consistent with a SFR that peaks at z ~ 1-2. 

To complement our treatment, we follow the approach 
developed in White & Ghosh (1998), Ghosh & White 
(2001) and Ptak et al. (2001): under the assumption that 
the X-ray luminosity of a galaxy results from the inte- 
grated luminosity of its populations of High and Low Mass 
X-ray Binaries (HMXB and LMXB), it is possible to de- 
rive the evolution of the X-ray luminosity by consider- 
ing the galaxy star formation history, the number ratio of 
HMXB and LMXB, and the timescales needed to switch 
on the X-ray emission. 

In order to derive a prediction for the evolution of an 
X-ray luminosity function, Ghosh & White (2001) and 
Ptak et al. (2001) considered different models of the CSFH 
in the redshift interval z = 0-8 representing the different 
determinations of the CSFH. Here we consider the two 
models discussed in Ptak et al. (2001): 

— a 'peak' model, which presents a rise of the SFR up 
z <~ 1.5 followed by an exponential decline, closely 
resembling the CSFH determination by Madau et al. 
(1998) (hereafter: 'peak-M' model); 
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Fig. 8. (In colour only in the electronic edition) X-ray counts derived from the IRAS luminosity function. The continuous and 
dashed curves, which converge at fluxes > 1CP 15 erg s" 1 cm -2 , show the number counts from the integration of the T03 XLF 
(upper, blue curves: integration and evolution performed up to z = 2; lower, green curves: integration up to z — 2 with evolution 
stopped at z = 1; continuous curves: 771 = 3.4; dashed curves: 771 = 2.7). Other symbols as in Fig. (7). Left panel: source counts 
from the T03 XLF with pure density evolution (77a = 3.4). Right panel: same XLF, but with pure luminosity evolutions (771 = 3.4, 
continuous curves and 771 = 2.7, dashed curves). 



— a 'peak+ Gaussian' model, which adds a Gaussian com- 
ponent at z = 1.7 to the 'peak' model, in order to ac- 
count for high redshift data, namely the sub- mm obser- 
vations by Blain et al. ( 1999a, b) (hereafter: 'peak-G' 
model). 

We derive X-ray number counts by considering the T03 
XLF and applying a pure luminosity evolution from the 
models in Ptak et al. (2001). The results are shown in 
Fig. (9). We notice that our predicted counts are in agree- 
ment with those derived in Ptak et al. (2001) by converting 
number counts in the Hubble Deep Field from the B band 
to the X-rays, for fluxes in the 0.5-2 keV band < 3 • 10~ 17 
erg s _1 ; altogether, the larger flux range (< 10~ 14 erg 
s _1 ) here explored the allows a direct check against the 
observed source counts. 

We find that the model of the XLF evolution derived 
from the 'peak-G' model of the CSFH, applied to the T03 
XLF, does not agree with the observed X-ray Log N- 
Log S. On the other hand, the 'peak-M' model is a good 
description of the observed source counts and closely re- 
sembles the source counts we have obtained by integrating 
the T03 XLF up to z = 2 with i]\ = 2.7 and stopping the 
evolution at z = 1. Thus it seems that the CSFH de- 
rived from sub-mm observations cannot be easily recon- 
ciled with the X-ray number counts. It might be, for ex- 
ample, that the sub-mm population of galaxies is obscured 
in the X-rays, or that a significant fraction of these galax- 
ies hosts an AGN, hence the true SFR at high redshift 



would be smaller. This latter hypothesis is also reinforced 
by the fact that the 'peak-G model' of the CSFH leads to 
a predicted value of the flux of the cosmic IR background 
~ 3 times larger than the current upper limit, as shown 
below. However, further consideration of this issue would 
be beyond the scope of this paper. 

Additional constraints on the evolution properties here 
examined come from the consideration that total flux of 
the galaxies should not exceed the cosmic background 
levels. The most stringent one comes from the Cosmic 
Infrared Background (CIB; see Hauser & Dwek 2001 for a 
review) since it is believed that star forming galaxies make 
the bulk of the CIB, while they only make around 50% of 
the radio background (Haarsma & Partridge 1998), or a 
few percent of the X-ray background (RCS03). Recent re- 
sults have shown that AGN contribute ~ 2% of the CIB 
at 60 n (Silva et al. 2004). 

Due to the extreme difficulty of the 60/i CIB level mea- 
surement, its actual flux value is still debated. A tenta- 
tive detection (Finkbeiner et al. 2000) fixed the flux value 
at vl v — 28 ± 7 nW m~ 2 sr _1 , while two upper lim- 
its have been set at 12 and 10 nW m~ 2 sr (Kashlinsky 
& Odenwald 2000 and Dwek 2001, respectively). Miville- 
Deschenes et al. (2002) also fixed an upper limit at 13.5 
nW m~ 2 sr -1 , with a tentative detection at 9 nW m~ 2 
sr -1 . While we refer the reader to Hauser & Dwek (2001) 
for a complete discussion of the CIB determination, in the 
following we will assume as the CIB flux 10 nW m -2 sr -1 , 
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Fig. 9. (In colour only in the electronic edition) X-ray counts 
derived from the IRAS luminosity function (Takeuchi et al. 
2003, 2004) with pure luminosity evolution according to Ptak 
et al. (2001). The solid curves show the integrated counts as- 
suming the 'peak-M' model for the cosmic star formation his- 
tory, with integration and evolution performed from z — to 
2 (lower, blue curve) or 4 (upper, pink curve). Dashed curves: 
same as above, for the 'peak-G' model. Other symbols as in 
Fig. (7). 

corresponding to the most stringent upper limit. Thus, the 
following estimates of the contribution to the CIB should 
be regarded as lower limits. 

The cosmic background flux from the integration of 
known sources may be defined as 

F= / dLogS— — S (8) 
J_ OQ dLogS 

of the differential Log iV-Log S, which in turn is related 
to the LF by 

dN /" Zmax dV" 
-P- = / dz V (LogL(S, z)) K{z) — (9) 
dLogS J Zmin dz 

where K(z) is the K-correction, here calculated after the 
spectral model of M82 from the GRASIL code (Silva et 
al. 1998). 

Taking the T03 LF as a reference and assuming lumi- 
nosity evolution up to z = 4, we found the contribution 
levels to the 60/i CIB to be in the range 60%-90%, depend- 
ing on the details of evolution: the lower figure was ob- 
tained by considering rj\ — 2.7 and evolution up to z — 2, 
the upper one by considering 771 = 3.4 and evolution up 
to z = 4. Had we considered density evolution, the con- 
tribution range would have been 60%-110% (lower figure: 
?7d = 2.7, z cv = 2; upper figure: r/d = 3.4, z cv = 4), thus 
ruling out density evolution beyond z cv ~ 2. 



The main conclusion that may be drawn, is that the 
galaxy LF may not probably sustain an evolution as strong 
as r/d — 3.4 or 771 = 3.4 beyond z ~ 2. Milder evolution like 
771 = 2.7 might be sustained up to larger redshifts, with 
the caveat that future analysis may lower the flux level of 
the 60fi CIB. 

In a complementary way, by using the 'peak-M' and 
'peak-G' models introduced above as prescriptions for the 
evolution of the infrared luminosity of galaxies, and consis- 
tently with the results from the Log A^-Log S analysis, we 
find that the the contribution to the CIB by sources dis- 
tributed according to the T03 LF with the 'peak-M' model 
for evolution is about 60%, while the 'peak-G' model is 
ruled out as the predicted flux from integrated sources is 
~ 290% of the CIB (due to the bell-like shape of these 
models, these percentages do not depend significantly on 
the maximum redshifts of integration and evolution). 

6. X-ray number counts and deep radio surveys 

Although not directly related to the determination of the 
XLF properties discussed in the preceding Sections, it is 
of interest to consider the X-ray number counts that can 
be derived from the radio sub-mJy population, which is 
commonly believed to include high redshift star forming 
galaxies as a major component. A prediction of the X- 
ray number counts based on the radio Log A-Log S for 
the sub-mJy population was already derived in RCS03. 
However, it was not possible in RCS03 to compare the 
predicted counts with observed ones, since the latter were 
not available at that time. 

The radio Log AT-Log S converted to the X-rays are 
shown in Fig. (10) as the short- (Richards 2000) and 
long-dashed (Fomalont et al. 1991) lines. The agreement 
with the observed counts is very good, with the Richards 
Log A^-Log S closely following the counts selected with 
the X-ray/optical flux ratio, and with the Fomalont et 
al. Log A^-Log S resembling the counts selected with 
the Bayesian criterion by Norman et al. (2004). Since 
Fomalont et al. (1991) checked that the source density at 
4 /ijy (5 GHz) derived from fluctuation analysis was con- 
sistent with a simple extrapolation of the Log AT-Log S, 
we extended the corresponding X-ray counts in down to 
the flux level of the radio background fluctuations. This 
simple extrapolation suggests that galaxies may become 
more numerous than AGN at some flux below ~ 10~ 17 
erg s _1 cm~ 2 . 

From Fig. (10) it is seen that X-ray sources with fluxes 
below ~ 10 -15 erg s _1 cm -2 should have a radio flux 
below ~ 1 mjy This may be checked as follows: we se- 
lected, from the sources with Log(Fx- ra y/-Fopt) < — 1 pre- 
viously discussed, the subsample lying in the CDFN area 
(97 sources); this subsample was cross-correlated with the 
1.4GHz catalogue of the same area (Richards 2000). We 
found that all the X-ray sources with a radio counterpart 
have a 1.4GHz flux below 2 mjy, with most of them around 
0.1 mjy. It has been proposed that a major fraction of the 
radio sources building the 'sub-mJy population', which 
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Fig. 10. (In colour only in the electronic edition) X-ray counts 
derived from deep radio Log TV-Log S in the 0.5-2.0 keV band 
(see RCS03). The short-dashed (blue) lines represent the 1.4 
GHz Log JV-Log S (Richards 2000). The long-dashed (red) 
lines represent the Log N-Log S by Fomalont et al. (1991). 
The vertical dotted lines show the limiting sensitivities for the 
radio surveys. For easier reading, the horizontal scale is shown 
both in radio (upper) and X-ray (lower) fluxes, assuming the 
conversion factors of RCS03. Other symbols as in Fig. (7). 

accounts for the majority of the radio number counts be- 
low ~ 0.5 mJy (Windhorst et al. 1985) and about half of 
the radio cosmic background at 1.4GHz (Haarsma et al. 
2000), are star forming galaxies at redshifts 0.5 < z < 1 
(Windhorst et al. 1990). Thus our findings contribute to 
reinforce this view. 

7. Conclusions 

By adopting well known relationships linking the IR, ra- 
dio, optical and X-ray emissions of normal spiral and 
starburst galaxies we have transformed available LFs of 
these objects from samples selected at IR, radio and opti- 
cal wavelengths into corresponding X-ray luminosity func- 
tions (XLF). From the comparison of these XLFs among 
themselves and with the X-ray selected XLFs of Norman 
et al. (2004) it is concluded that: 

— A clear prediction for a local (z = 0) XLF emerges 
from the comparison of the infrared, radio and optical 
LFs (Fig. 3): the derived XLFs agree within a factor 
2 in the luminosity interval 10 40 -10 41 erg s _1 encom- 
passing the knee region after which all XLFs steepen 
toward higher luminosities; although departures at 
lower and higher luminosities are present, the average 



local X-ray luminosity density (~ (3 • 10 37 ± 30%) erg 
s" 1 Mpc^ 1 ) appears to be well defined (Table 2). 

— There is a remarkably good agreement between the two 
best defined XLFs from the IRAS sample (T03) and 
W03, both based on the largest galaxy samples with ~ 
10 4 objects each and on very different selection criteria: 
the W03 XLFs are all well reproduced by the local 
T03 XLF once a luminosity evolution with coefficient 
rji ~ 3 is applied. Thus the T03 XLF may be considered 
representative of the galaxies XLF up to z ~ 1 . 

— None of the XLFs derived from the infrared and radio 
local samples (z < 0.3) can match the higher luminosi- 
ties (Lx > 3-10 41 , Figs. 4a, c) of the observed (Norman 
et al. 2004) XLFs with their quoted (Table 1) den- 
sity evolution parameters (%)■ However, a very good 
match is achieved, particularly in the high redshift bin, 
by adopting pure luminosity evolutions with param- 
eters rji approximately equal to 77^ of the respective 
samples (Figs. 4b, d). Among the XLFs from the radio, 
IR and optical samples selected over a wider redshift 
range (P04, C89/H00, W03) only those from the blue 
sample (W03) match the observed XLFs in both red- 
shift bins and in the luminosity intervals in which they 
are defined. 

By combining the above findings we conclude that 
the local XLF derived from IRAS T03 sample provides 
a consistent description of the local XLF of star form- 
ing galaxies and that pure luminosity evolution of the 
form (1 + z) m , with r/\ ~ 3 (consistently with i]\ ~ 2.7 
as found by Norman et al. 2004, and within errors of the 
ijd = 3.4 ± 0.7 coefficient given for density evolution by 
T03), is favoured over density evolution providing a good 
representation of available data as derived from radio, op- 
tical and X-ray samples. We note, however, that the excess 
of the observed XLF at higher luminosities in the low red- 
shift bin may be caused by a slight pollution by AGN in 
the Norman et al. (2004) sample. 

By comparing the number counts from our adopted 
XLF from the T03 sample with those derived from ob- 
served X-ray samples of star forming galaxies, and with 
the total X-ray number counts, we reach the following 
conclusions: 

— The integral number counts from samples of star form- 
ing galaxies show very similar slopes but their normali- 
sations are dispersed within a factor ~ 2. The observed 
fraction of star forming galaxies at the flux limit of the 
Chandra Deep Field surveys (~ 3-10 -17 erg s _1 cm -2 ) 
is (30 ± 12)%. 

— The counts obtained by adopting different values of the 
luminosity evolution parameter 771 and redshift cut-offs 
in the evolution indicate that a good representation of 
the observed data is obtained with 771 ~ 2.7 and a cut- 
off of the evolution in the redshift interval 1-2, which 
is believed to host the peak of the star formation rate. 
Much stronger evolution, say rji — 3.4 as adopted in 
our analysis, is not consistent with the observed counts 
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unless one stops the evolution of the XLF at a redshift 
z ~ 1. 

— Additional, supportive evidence of the above results is 
provided by the very good agreement of the X-ray inte- 
gral counts derived from the radio Log iV-Log S of the 
sub-mjy population, believed to be mainly represen- 
tative of the star forming galaxies, with the observed 
counts and those from the integration of the T03 XLF 
(compare Fig. 7 and Fig. 10). 

We have complemented the analysis of the Log TV- 
Log S by adopting our proposed XLF in the framework 
of the cosmic star formation history (CSFH) as parame- 
terised in the models of Ptak et al. (2001), finding that the 
observed counts are only consistent with model 'peak-M' 
of these authors, i.e. with a model of the CSFH resem- 
bling the data from Madau et al. (1998). We note that the 
CSFH model suggested by sub-mm observations (Blain et 
al. 1999a) appears to be ruled out both by the observed 
X-ray number counts and by the flux level of the Cosmic 
Infrared Background (CIB). 

As a further check on our results, we have estimated 
the contribution of the galaxies in the T03 sample to the 
CIB. Pure luminosity evolution with rji <~ 2.7 is consistent 
with the current CIB upper limit; moreover, the X-ray 
counts predicted by the 'peak-M' model are almost coin- 
cident with those obtained from the integration the T03 
XLF with r]\ = 2.7. We have also found that the evolu- 
tion implied by the 'peak-G' model on the T03 LF would 
exceed the CIB limit by a factor ~ 3. 

It may be predicted that 6 • 10 3 ~10 4 sources deg -2 
would be present with fluxes > 10~ 17 erg s~ 4 cm~ 2 . This 
figure may be compared with the AGN number density 
between 4 • 10 3 and 10 4 deg~ 2 at fluxes below <~ 10~ 16 erg 
s~ 4 cm~ 2 , as predicted by AGN models (Comastri et al. 
1995; Ueda et al. 2003). It follows that, if we conservatively 
assume the upper estimate, the X-ray number counts of 
normal galaxies should overcome the counts from AGN at 
fluxes < 10~ 17 erg s~ 4 cm~ 2 (6- 10~ 18 erg s~ 4 cm~ 2 in the 
case of the W03 XLF or the T03 XLF with the 'peak-M' 
evolution) . 

The Chandra Deep Field surveys have provided for the 
first time a sizable number (100-200, depending on selec- 
tion criteria) of X-ray detections of high redshift galax- 
ies. However, further studies on galaxy evolution would 
greatly benefit from enlarged and/or deeper samples. If a 
0.5-2 keV flux of 10~ 17 erg s~ 4 cm~ 2 (<~ 3 times fainter 
than the current flux limit of the Chandra Deep Field 
North) could be reached, it would be possible to detect 
galaxies with Lx = 10 41 erg s _1 placed at z = 1.3, or 
with Lx — 3 • 10 41 erg s _1 at z = 2. Our models show that 
the density of galaxies within this luminosity and redshift 
ranges should be <~ 4000 deg -2 ; that is considering over 
an area of, say, 100 arcmin 2 , ~ 100 could be detected. In 
principle it would then be possible to derive the high lu- 
minosity tail of the galaxies XLF in the z = 1.3-2 redshift 
bin, thus placing a significant constraint on the evolution 
of the SFR. 



Although these studies are feasible in principle with 
Chandra, they would require a huge investment of time. 
From our models it is also apparent that more stringent 
constraints could be achieved by pushing the detection 
limit toward 10~ 18 , so that XEUS would be more suited 
for this purpose, if the goal of a PSF with diameter < 2" 
can be achieved (given the predicted number of galaxies 
around 10~ 18 erg s" 1 cm -2 , a significantly larger PSF 
would lead to source confusion). 
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